Robustness of Majorana fermions in proximity induced superconductors. 
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In 2D chiral p-wave superconductors, the zero-energy Majorana fermion excitations trapped at vortex cores 
are protected from the thermal effects by the mini-gap, A'^/ef (A: bulk gap, €f' fermi energy), which is the 
excitation gap to the higher-energy bound states in the vortex cores. Robustness to thermal effects is guaranteed 
only when T <C A^/ef ~ 0.1 mK, which is a very severe experimental constraint. Here we show that 
when s-wave superconductivity is proximity-induced on the surface of a topological insulator or a spin-orbit- 
coupled semiconductor, as has been recently suggested, the mini-gaps of the resultant non-Abelian states can 
be orders of magnitude larger than in a chiral p-wave superconductor. Specifically, for interfaces with sufficient 
barrier transparencies, the mini-gap can be as high as ~ A 3> /ep, where A is the bulk gap of the s-wave 
superconductor. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45.+C 



I. INTRODUCTION: 

Topological quantum computation (TQC) requires the ex- 
istence of a 2D topologically ordered state whose lowest- 
energy excitations follow non-Abelian statistics'. If the ap- 
propriate many-body ground state wavefunction - e.g., Pfaf- 
fian states in fractional quantum Hall systemsi*^ and chiral 
p-wave (px + ipy) superconductor/superfluid^'"' - is a linear 
combination of states from a degenerate subspace, then a pair- 
wise exchange of the particle coordinates can unitarily ro- 
tate the wavefunction in the degenerate subspace. This exact 
non-Abelian statistical property can be used to perform quan- 
tum gate operations, which are, in principle, fault-tolerant''^'^. 
More importantly, these non-Abelian particles, the Majorana 
fermions, being half-fermions, are new particles in nature dis- 
tinct from ordinary Dirac fermions, which are of obvious in- 
trinsic fundamental interest^,. 

In practice, a key requirement for TQC is that the degen- 
erate ground state subspace must be separated from the other 
excited states by a non-zero energy gap, so that thermal effects 
cannot hybridize the topological quasiparticle states with the 
other higher-energy, non-topological, states in the system'. In 
2D Px+ipy superconductors (SC), such as possibly SrRuOi, 
where the zero-energy Majorana fermion excitations trapped 
in the vortex cores are the topological quasiparticle states, this 
gap is given by the so-called mini-gap, ^ (5o ^ A^/ej?, where 
A is the bulk superconducting gap and e^? is the fermi energy^. 
Since (5o < 0.1 mK is a very small energy scale for typical 
p-wave superconductors, the requirement T ^ do constitutes 
the real bottle-neck for TQC, even if the best possible 2D p^ + 
ipy superconductor-based platform were realized in the labo- 
ratory. Here we show that, in a class of newly-proposed TQC 
platforms, involving Majorana fermions in multilayer struc- 
tures where s-wave superconductivity is proximity-induced 
on a host topological insulator (TI)2iil or a spin-orbit-coupled 
semiconductor'"*, the mini-gap can be enhanced by several 
orders of magnitude. Given that a strong proximity effect 
in such superconductor-semiconductor structures has already 
been experimentally demonstrated'^ '^, it is realistic to de- 



crease T to satisfy T <C So, since Sq can be made as high 
as 1 K (i.e. a fraction of order unity of A), which is the bulk 
gap in the s-wave superconductor. 



To derive these results, we explicitly analyze the micro- 
scopic model of the proximity effect between a TI surface 
and an s-wave superconductor by applying the conventional 
tunneling formalism''^. We find that, in addition to the super- 
conducting gap A, the interface transparency (denoted by A 
below) given by the inter-layer tunneling amplitude controls 
the strength of the proximity effect on the TI surface. Our 
main result is that for barriers with transparency satisfying 
ef ^ A [/, A, where U is the fermi level on the TI surface, 
the excitation gap above the non-Abelian quasiparticle states 
can be ~' A;» A^/ej?. We show this by applying our central 
result, Eq. (13), on the excitation gaps in the two most impor- 
tant structures on a TI surface, a line junction (Eq. (16)) and a 
vortex (Eq. (19)), which contain Majorana modes. Note that, 
as discussed earlier'', the Majorana modes in a line junction 
and a vortex are the two most essential elements of a putative 
TQC architecture on the TI surface. The dramatic increase 
of the excitation gap on the TI surface greatly enhances the 
robustness of the topological quasiparticles to thermal deco- 
herence effects, which may bring non-Abelian statistics and 
TQC to the realm of realistic, achievable, temperature regimes 
in the laboratory. 



The paper is organized as follows. In section II, we describe 
the microscopic model we consider for the proximity effect at 
a TI-SC interface. In section III, we derive the proximity in- 
duced effective pair potential with renormalized parameters. 
In sections IV and V we described the line-junction geometry 
that may be used to manipulate superconducting quasiparti- 
cles such as Majorana fermions at a TI-SC interface. In sec- 
tion VI, we show that the mini-gap of a vortex formed at a 
TI-SC interface can be orders of magnitude larger than in a 
vortex in an intrinsic (not proximity induced) superconductor 
such as SrRuO^. 



II. MICROSCOPIC MODEL FOR PROXIMITY EFFECT 

An interesting property of superconductors is that they can 
induce superconductivity in a normal metal in contact with 
the superconductor.'^ This is referred to as the superconduct- 
ing proximity effect. As we will show in the rest of the text, 
the superconducting proximity effect allows a greater degree 
of control of the superconducting quasiparticle spectrum, than 
is possible by simply modifying the superconductor, where 
the superconductivity is intrinsically derived from the quasi- 
particle spectrum itself. The superconducting properties of a 
normal-superconductor (NS) interface can be described by the 
self-consistent Bogoliubov-de Gennes (BdG) equations at the 
interface. 

The BdG equations at an NS interface can be writ- 
ten in terms of a Nambu spinor wave-function, 5'(r) = 
{uf{r),ui{r),vi{r), —v^{r))^ which is finite on both the su- 
perconductor and normal side of the interface. In order to 
write these equations more compactly it is convenient to in- 
troduce the Nambu matrices Tx,y,z, which are identical to the 
Pauli spin matrices (Jx,y,z, except that they operate on the 
(m, v) part of the spinor ^'(r). Thus the spinor '^{r) exists 
in the tensor product space aa® Tp. The BdG equations for 
the quasiparticle wave-functions can be written in terms of 
the a ® T matrices in the form Hsdc^i''') = E'^{r), where 
HbcLG is a 4 X 4 BdG Hamiltonian. 

The normal-superconductor interface can be considered to 
be a planar geometry (Fig. with the coordinates r = 
(r, z) where r and z are the in-plane (r = (x, y)) and out- 
of-plane coordinates for the interface (the N-S interface is at 
z — 0). In this paper, we will restrict ourselves to the case 
where the superconductor in our system (SC) is a conventional 
s—wave superconductor such as Al, and the normal part(N) of 
our system is the metallic surface state of a topological insula- 
tor (TI) such as i?i2re3.'" The BdG Hamiltonian describ- 
ing the s—wave superconductor is of the form 

^sc- ^^-EF T,+A,(r)r, (1) 

\ 2m* / 

where the m* is the effective mass of the superconductor, ep 
is the fermi energy of the superconductor and l^sif) is the 
pairing potential in the supercondutor. The pairing potential 
As(r) inside the s—wave superconductor is generated by a 
pairing interaction V{r) that can be taken to be point-like.'^ 
In the self-consistent BdG equations, the pairing potential sat- 
isfies the constraint 

A,(r) - V{r) <,t(r)«„4(^)(l - 2/(i?„)) (2) 

n 

where (u„,w„) are components of the 4-spinor eigenvec- 
tors 5'„(r) which satisfy the BdG equations Hsc'^,i{r) = 
-£'„^'„(r). Here f{En) is the fermi-function. 

The normal part of the TI-SC system, which in our case 
consists of the metallic surface bandiS^i^ of a topological in- 
sulator (TI)(in anticipation of the proximity induced super- 
conductivity) can be described by a BdG Hamiltonian 

HtI^IiVCT ■Vr-U]lg)Tz, (3) 




FIG. 1 : (Color online) Proximity induced pairing on the TI surface. 
The (red) region on the left is the topological insulator (TI), and the 
(blue) region on the right is an s-wave superconductor. 'Integrating 
out' the superconducting degrees of freedom produces the self en- 
ergy E on the TI surface, where E is given by the tunneling Hamilto- 
nian T and the Green function Gqq of the superconductor (see text 
for details). 



where U — ep — J dz\<f>{z)\'^Vgate{z) is the fermi level of the 
TI surface where Vgate (z) is the gate potential and (f>{z) is the 
z-dependent electron wavefunction (with momenta close to 
the Dirac point) of the TI surface states. Here v is the effective 
electron velocity on the TI surface. Note that apart from the 
factor Tz, this is the Dirac-Hamiltonian describing the surface 
state of the topological insulator in the normal state. The fac- 
tor Tz is present to account for the independent electron-part 
(represented by u in the spinor ^P) and hole-part (represented 
by V in the spinor 'I'). 

Since there is no tunneling so far between the TI and SC, 
the TI may be described as being normal with no supercon- 
ductivity in the TI. The introduction of a tunneling term T 
which transfers electrons between the SC and TI leads to a 
finite value for the order parameter {ilj„{r)tjja-' {r')) on the TI 
surface.'^ However it is crucial to note that despite the exis- 
tence of an order parameter in the normal part of the system 
{TI), the microscopic pairing potential At/(t) = 0. This 
follows from the self-consistency condition in the TI, 

Ariir) = VTi{r){^f{r)iJj^{r)) = (4) 

since the attractive pairing interaction in the TI van- 
ishes (Vriir) = 0). This is consistent with the de 
Gennes boundary conditions at the interface which re- 
quires A{r) / N {r)V (r) to be continuous across the TI-SC 
interface.'^ Here A(r), V{r), N{r) stand for the microscopic 
pairing potential, pairing interaction and density of states at 
the fermi level on both sides of the interface. Note that since 
Ati and Vri are both zero on the TI side of the interface, the 
ratio can be finite which allows tt^^t — = ,f^.^, . However, 
even though there is no microscopic pair potential on the TI 
side, the superconducting proximity effect induces an effective 
superconducting pair potential. In the next few paragraphs, we 
show explicitly how the superconductor can be integrated out 
to yield an effective Hamiltonian on the TI surface which has 
an effective pairing potential. 
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III. EFFECTIVE PAIRING POTENTIAL ON THE TI 
SURFACE 

The complete BdG Hamiltonian descrbing the TI-SC inter- 
face including the tunneling term'*^ (Fig. ([T]i) is defined by 
the Hamiltonian: i?totai = -f^Ti + ^sc + T + T^. Here, 
Hti and Hsc are the Hamiltonians describing the TI surface 
and the s-wave superconductor, respectively. T describes the 
tunneling from the TI surface to the superconductor and 
describes the tunneling in the opposite direction. The exci- 
tation spectrum of the interface can be determined from the 
Bogoliubov-de Gennes (BdG) equation 



(iJtotai - E)^ir) = 0, 



(5) 



where 5'(r) is the Nambu spinor , and iJtotai is written as a 
4x4 matrix in the Nambu basis. 

The tunneling Hamiltonian Ht = T + T'^ coupling the 
2D TI surface states with the superconductor can be explicitly 
written in the Nambu space as 

T{r;r'z') = t,J d^kdk.x{^';kk,){x{k)\r\(t>}e'''<^-'^'\ 

(6) 

Here the momenta are measured relative to the Dirac cone mo- 
mentum M and the tunneling matrix element in the integrand 
can be approximately written as^ 

(x(fc)|T|0) = — [(j){z)dzx{z; k, k^)~x{z; k, k^)d:,(l){z)]\^=f), 
m 

where, fc, kz) is the single-particle eigenfunction in the 
superconductor and 0(z) is as defined before. 

In order to solve the BdG equation at the TI-SC interface, 
we decompose the wave-function as = ip^i + i/isc- De- 
composing the BdG equation (Eq. (|5])) we obtain 



(i/Ti - s)^Ti + rVsc - 



(7) 
(8) 



Substituting the wave-function i/^sc from Eq. ([H]) in Eq. (|2| we 
get the effective BdG equation on the TI surface, 

{Hti + S(r/; uj) - w)i/'ti = 0. (9) 

Here the self-energy S on the TI surface (Fig ([T]i) is given by 

E(rr',w) = drirfr2r^(r,ri)4c(ri,r2;c^)r(r2,r'), 

(10) 

where Gg^(ri, r2; lo) — (i?sc — t^^)^^ is the Green function 
matrix in the superconductor. 

The effective superconducting pairing potential induced at 
the surface of a TI appears in the form of the anomalous part 
of the self-energy, S, which can be written as 

T,A{rr\uj) = - J dridr2T^r,ri)F^°(}{ri,r2;uj)T{r2,r'), 

(11) 



where F^°}{ri,r2;uj) = J dte"^* {Ti}{ri,0)^p{r2,t)) is the 
anomalous part of the Green function in the superconductor 
that represents the superconducting order parameter. Since 
the anomalous part of the self-energy is proportional to the 
anomalous part of the Green function, the complex phase of 
Eyi must also equal the complex phase of the superconductor. 

The Green function for the superconductor can be written 
in terms of the normal state eigenbasis of the superconductor 
as, 

4°(] (rr'; c.) = ^ Xn{r}Xn {r'K{e,, - 8^^, + At, - 

n 

(12) 

where the normal state eigenfunctions Xn are taken to be real 
and spin independent. Using Eq. ^ and Eq. ( fT2] i in Eq. ( fTOl i 
and then Fourier transforming to the momentum space the 
self-energy on the TI surface takes the form, 

S(rr';a;) = -^X«(n)x«(r2)T'''(r, r^)T{r2,r') 

n 

wTo + (e„ - eF)Tz + Atx 



(13) 



This equation is more conveniently written in terms of the tun- 
neling density of states A on the TI as 



I]{rr' ;ijj) = —J de 



2^ A2 —Hrr';e) (14) 



where the tunneling density of states A from the TI to the SC 
relative to the fermi level is given by 

A(rr';e) = ^ J(e-e„+£)x„(ri)x„(r2)rt(r, ri)r(r2, r'). 

n 

(15) 

Note that the matrix A is the operator that can be used to cal- 
culate the time a state (l){x) on the TI would take to diffuse into 
the SC layer when in the normal state through the expression 
T-i = / drdrXr)^*(r')A(rr';e)- 

The scale of the spatial and energy dependence of the ma- 
trix A{rr';e) can be understood by observing that it is de- 
pendent only on the wave-functions Xn of the superconduc- 
tor in the normal state at the energy e. Defining the Fourier 
transform of A(rr'; e) as A(r, k;e) — J dr'A{r + r'/2, r — 
r'/2; e)e^''"' , similar to the quasi-classical approximation, 
we obtain a function in position r and wave-vector k space. 
Given that the original definition of A involved only the nor- 
mal state bandstructure of the superconductor, A{rk; e) varies 
on a scale k ^ kp and € ^ ep- Since we are interested in slow 
variations ^ kp and energies J7, A <C ef, we can ignore 
the k and e dependence of A and take it to be a function of 
only r. Thus we will assume A(rfc; e) w A(r) — A{rkp; ep). 
Fourier transforming back to real space leads to the real space 
relation A(rr'; e) ~ X{r)S{r — r') which is local in r, r' and 
independent of e. Within this approximation we find. 



{-UJTQ + ATj;) 

x/|AP-a;2 ■ 



(16) 
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Using Hti from Eq. (2) and the local (fc-independent) self 
energy from Eq. ( fT6] l. we can now straightforwardly rewrite 
Eq. ^ as an effective BdG equation for the TI surface: 

[ij{uj)tcr ■ Vt^ - U{uj)t^ + A(w)t^ - wJV'Ti = 0, (17) 

where = Z{lu)v, U{uj) = Z{uj)U and A(w) = 

AA/(VA2 -a;2 + A). Here, the factor = 
VA2 - w2/(VA2 - a;2 + ^j-^)^ ^nd A(w) are the 
renormalized velocity, fermi level, and superconducting gap 
on the TI surface, respectively. Below we will be interested 
only in the low-energy states with energies E A (typically 
0.1 A < E < 0.5 A). In this case, we can approximate the 
frequency-dependent parameters in Eq. ([TtI i with their values 
at w = 0: 

1+^ -"-'A -"-'A 

(18) 

The renormalization of the parameters described in Eq. ( fTSl l 
gives the central results of this paper which can be understood 
as arising from the virtual propagation of the electron in the 
superconductor. This is consistent with the estimate of the 
fraction of time spent in the TI, which is tri / {tri + tsc) ^ 
(1 + A/A)~^ = Z{{)), since we can estimate life-time of an 
electron in the TI as txi ~ A^^ and the virtual excitation of 
an electron in the SC as tsc ^ A~^. Thus, in addition to 
inducing a pairing potential A(aj) on the TI surface, the prox- 
imity effect also renormalizes the velocity on the surface of 
the TI to w — > v{uj) = Z{uj)v and the background potential 
to [/ — > U{uj) ~ Z{lu)U. One can notice that the TI electron 
density n is proportional to (U/v)'^ and remains constant as 
expected. In the weak tunneling regime \ ^ ep, our tunnel- 
ing matrix approach to the proximity effect is consistent with 
previous calculations^" for the superconductor-semiconductor 
system using the Blonder- Tinkham-Klapwijk model^'. Since 
the parameter A is determined by purely electronic energy 
scales, A can be larger than A and the retardation effects dis- 
cussed above lead to substantial renormalizations of the orig- 
inal parameters. 

The expression for the proximity-induced self-energy, 
Eq. [18] was derived for a single superconducting island in 
contact with a TI surface. The proximity induced pairing was 
then localized to the region where A(r) was non-zero. In the 
case of different superconducting islands covering different 
parts of the TI surface, one must add the self-energy contribu- 
tion from each of the separate islands. The pair-potential A 
is then taken to be a constant for each island, but varies be- 
tween different islands. In particular, a region of the TI which 
is in contact with a superconducting island (through a non- 
zero A(r)) inherits a superconducting pairing potential with a 
phase which has the same value as the island it is in contact 
with. Futhermore the effective pair potential, vanishes (since 
A(x) = 0) if the TI surface is not locally in contact with the 
superconductor (for example the TI regions in between two 
superconducting islands has vanishing pair potential). 

Below we will apply the formulae contained in Eq. JTSl l to 
estimate the excitation gaps above the Majorana topological 
excitations that have been discussed for the TI surface -. 



rV. MAJORANA SYSTEM ON TI SURFACE 

The proximity effect discussed in the previous sections al- 
lows one to induce superconductivity on a semiconducting 
substrate such as a TI surface state. This leads to the pos- 
sibility of manipulating the position and energies of super- 
conducting quasiparticles on such a surface state. One of the 
most interesting scenarios that this kind of manipulation of 
superconducting quasiparticles might lead to is the possibil- 
ity of creating and braiding Majorana fermions at such TI/SC 
interfaces. 

The trijunction geometry shown in Fig. |2] provides a way 
to create and braid a Majorana fermions on the TI surface: 
two tri-junctions (A and B in Fig. |2]i of superconducting lay- 
ers with distinct local phases separated by a line junction of 
length L . The superconducting islands can be considered to 
have distinct phases for barrier transparencies much less than 
unity corresponding to A <C ep- The line-junction separating 
the superconducting islands consists of a region of width W 
which is not in direct contact with either of the superconduct- 
ing islands. Therefore the effective pairing potential vanishes 
in this region and in principle the line-junction can support 
multiple low energy transverse modes with energy much less 
than A. This can be avoided by taking the junction width to 
satisfy W ~ v/ A = v/X such that the confinement energy 
v/W creates a sufficiently large energy spacing between the 
various transverse modes. For the transparencies assumed in 
the paper A ^ A <C ef, the width W , which is required to 
be smaller than the effective coherence length is much smaller 
than the coherence length of the superconductor in the clean 
limit but much larger than the fermi wave-length of the super- 
conductor Ideally, one should adjust the barrier transparency 
to set A to the lowest value such that the line-junctions width 
W is within the resolution of the fabrication technology. 

Such a geometry can be used to trap zero-energy Majorana 
fermions on the TI surface at discrete vortices formed by the 
phase configurations of the superconducting islands either at 
A or B. By varying the phase differences between the islands 
it is also possible to move the discrete vortex from A to B 
and thus transport the bound Majorana state. The phase dif- 
ferences between neighboring superconducting islands is con- 
trolled by connecting the islands by superconducting loops 
with fluxes threaded through them. 

For barrier transparencies between the TI and the SC, such 
that A ^ A, but still much less than unity (which implies 
A <C ep), the Josephson current between the islands mediated 
by the single channel TI surface state splits up into a small cur- 
rent density in the bulk-like superconducting loops. This adds 
only small gradients in phase to the gradient resulting from the 
vector potential of the magnetic flux. Thus the difference in 
the phases on the islands are controlled by the applied fluxes 
according to the equation 

- 0fc = $jfe. (19) 

Therefore, for the application of the fluxes shown in Fig. |2] 
leads to the desired phase configuration on the islands. 

The tri-junction set-up discussed in the previous paragraphs 
operates on the principle that if one considers a single line- 
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FIG. 2: (Color online) A tri-junction-pair geometry of superconduct- 
ing islands deposited on the TI surface (top view) to confine and ma- 
nipulate Majorana fermions. The superconducting islands are con- 
nected strong superconducting loops enclosing fluxes ^n,m- The 
values of the fluxes $1^2 = $3,4 = f + 50, $2,3 = f — 50, 
$4,1 = — 3|- — S(f) which satisfy — control the phases 

of the superconducting island. For the given superconducting phase 
configuration and 5<p = j the structure contains a vortex with a 
trapped Majorana state only on tri-junction A. By changing Stf) to 
the discrete vortex together with the Majorana state is trans- 
ferred to the tri-junction B. The Majorana state is transported from 
A to B by the delocalized Majorana fermion state formed on the ID 
line junction (of length L) connecting A and B in the intermediate 
stage with S<l) — 0. 

junction separating 2 superconductors (such as the one join- 
ing A and B) with a phase difference of tt, there are exactly 
1 pair of energy and momentum Majorana states trapped 
between the islands. At finite but small momenta, one can use 
lowest order k.p perturbation theory to check that these modes 
disperse into a pair of linearly dispersing non-chiral Majo- 
rana modes with equal and opposite velocities. Off-setting the 
phase difference from tt by 6(f) leads to a mass gap in the lin- 
ear Dirac spectrum of this pair of Majorana modes. Thus the 
effective Hamiltonian describing such a line-junction is given 
by 

Hjunc = -1-^ TTT'^z^ + m(H)<^x (20) 

where a here is an effective pseudo-spin degree of freedom 
that keeps track of the mode index. The term proportional 
to m is an effective phase difference dependent mass term. 
Following the Jackiw-Rebbi index theorem^^ it is clear that a 
localized zero-energy solution results at any point where the 
sign of the mass mScj) changes sign. Such changes in signs can 
be engineered at the ends A and B of the line-junction joining 
A and B by appropriately tuning the fluxes. As it turns out, 
the appropriate mass changes are generated at the ends A and 
B whenever there are discrete vortices in the phases at A and 
B^ 

With the configuration of the phases on the superconducting 
islands as shown in Fig.|2] and for 6(f) = ^, the total phase 
change around the tri-junction A is 27r. Tri-junction A then 
acts as a discrete vortex and there is a localized zero-energy 
Majorana state confined to A. In this configuration, there is 



no vortex or zero-energy mode at B. It can be easily checked 
from Fig.|2]that the roles of A and B are reversed if 6(f> = 
now B contains a vortex and a localized Majorana mode while 
A is topologically trivial. In both cases, the spectrum of the 
line junction connecting A and B is gapped with the excitation 
gap controlled by d(p. To avoid hybridization of the localized 
states at A and B, the length L must exceed the size of the 
locahzed states themselves, 

L>^r^ i/A. (21) 

where is the decay length of the Majorana states on the TI 
surface. 

It is now clear that the Majorana states trapped at the dis- 
crete vortices can be braided by tuning the phase 6(fi through 
zero. For 6(f> = the phase change across the line junction is 
TT, and for the arrangement of the phases as shown in Fig. |2] 
there is a single zero energy extended Majorana mode on the 
line junction. When 5(t> is tuned from -I to to — -I, the Majo- 
rana mode shifts from A to the line junction and finally to B. 
For 6(f) = 0, the other low-energy delocalized modes on the 
line junction follow a dispersion given by^, 

u{q)^±qiAy{U^ + A^), (22) 

Below we will consider two types of excitation gaps which 
control the thermal robustness of the above Majorana sys- 
tem. First, in the line junction of length L the gap Eg w 
I A2/(C/2 + A^) that follows from Eq. ^ protects the de- 
localized zero-energy Majorana mode from thermal decoher- 
ence. Eg controls the thermal robustness of the Majorana 
fermions while they are braided in TQC. We show below by 
explicit analytic arguments that it is possible to make Eg ^ A 
by appropriately designing the TI-SC interface. The thermal 
robustness of the (stationary) topological qubits themselves, 
on the other hand, is determined by the energy gap {AE) 
above the zero-energy localized Majorana states within the 
discrete vortex cores. We will show by rigorous numerical 
calculations that even this scale AE ^ A, making the entire 
TQC architecture surprisingly robust to thermal decoherence 
effects. 



V. EXCITATION GAP IN LINE JUNCTION 

For a line junction of length L, the gap Eg is given by (see 
Eq. (|22li) Eg w |:AV(C/2 + A^). Now, for barriers with 
transparency such that A ^ U,A, we get A ^ f/ (Eq. ( fTSl l) 
and the factor multiplying v/L in Eg reduces to unity. Even 
\fU^\, which should be possible experimentally, this fac- 
tor is still of order unity. To maximize Eg, we need to take 
the minimum allowed value of the length of the line junction, 
Lm ^ v/ A (Eq. (I2TI) ). Therefore, the maximum E'g attainable 
on the TI surface is given by, 

Eg^ -^^A= (23) 

Lrn 1 + ^ 

which, in the case of transparency A ^ A, reduces to A itself. 



VI. EXCITATION GAP IN VORTEX: 

To determine the excitation gap Ai? within a vortex core 
numerically, we consider the BdG Hamiltonian on the surface 
of a TI sphere with a vortex and an anti-vortex at the poles^^, 

= [wi? • (a X p) - U]t, + A(r)T, (24) 

which can be written in angular coordinates as 

H = [--^I/-CT-C/]r2 + A(6'){cos(/)r^+sin0-ry}. (25) 
it 

Here R is the radius of the sphere, A(0) = 
Atanh{i?sin(?/^„} and is the size of the vortex 
core. In the above Hamiltonian, we have approximated 
discrete vortices by regular ones with continuously vary- 
ing phases. The resultant azimuthal symmetry allows 
us to decouple the equations into sets indexed by m 
with a basis of spinor spherical harmonics of the form 
{Yi^r,i+i,yi,m+2,Yi^„i,Yi^rn.+iY{0-,4>)- We cxpcct the 
minigap of such a continuous vortex to be qualitatively 
similar to the discrete vortex in Fig. (|2]l. We find that the 
m — —1 channel contains a pair of decaying and oscillating 
solutions which are spatially localized at the two poles. The 
corresponding eigen-energies exponentially decay to zero 
with the radius of the sphere, indicating that, in the limit 
when the vortices are far-separated, the eigen-energies are 
exactly zero. On the other hand, the spectrum of the other m 
channels qualitatively resemble the m = — 1 channel, with 
the important difference that the eigen-energy of the first pair 
of excited states does not vanish as the radius of the sphere 
increases. This eigen-energy gives us the excitation gap in 
the vortex core. 

Assuming the vortex core size to be equal to — = 
w/A — v/X, the numerical results for the mini-gap (Fig. 2) 
can be fit by the analytic form 

AE^O.mA^/VA^ + C/2. (26) 

The above result is interesting since it approaches the value 
AE ^ A ( avoiding the approximation AE ^ A in Eq. [18] 
yields AE sa 0.6A instead of 0.8A) for - A . Such a 
large mini-gap is made possible for proximity induced super- 
conductors because of the ability to tune the chemical poten- 
tial on the TI surface independent of the chemical potential of 
the superconductor This is unlike the case of more conven- 
tional s-wave and chiral p-wave superconductors where the 
chemical potential U ^ A and the minigap obeys the scaling 

A2 

AE - 0.83-^ (27) 

similar to the classic result of Caroli, de Gennes and 
Matricon.** 

In fact this is expected since the scaling of Caroli, de 
Gennes and Matricon can be established using fairly generic 
semi-classical arguments by considering the core of the vor- 
tex to be a normal region of radius ^ with no superconductivity 
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FIG. 3: (Color online) Numerical results for the vortex mini-gap 
AE plotted against the renormalized fermi level U on the TI surface 
{AE, U scaled by A). The solid (red) line gives the mini-gap when 
the vortex core size = ^ ~ v/A — v/X, as is appropriate in a 
regular vortex with a continuously varying phase. The dashed (black) 
line shows that the excitation gap is even larger when the vortex core 
size is smaller, as is expected in a discrete vortex (see Fig. 1). The 
inset shows the TI sphere with a vortex and an antivortex (with re- 
duced superconducting amplitudes at the vortex cores) situated at the 
north and the south poles. 



which can host quasiparticle states close to the fermi level ep. 
Considering a quasiparticle state at angular momentum that 
is confined in the vortex core, the state at the next allowed an- 
gular momentum would have a relative energy vp{\k\ — kp) — 
vp{y^kp + — kp). Here kt ^ is the transverse mo- 
mentum from the angular momentum. For kt <C kp, this 
leads to the energy splitting vpkf /kp ^ vp/kpS^"^ . For a vor- 
tex, the quasiparticles are confined within a decay length of 
^ = vp/ A. The leads to the relation 5q ^ A'^ /ep. Since this 
argument applies to any system with weak superconductivity, 
it also holds for the TI/SC system as verified by our numerics. 

Thus the staggering difference between the mini-gap esti- 
mate in a chiral p-wave superconductor (5o = lO^^meV ^ 0.1 
mK and the estimate of 1 K for the proximity induced system 
arises entirely from our ability to control the fermi energy ep 
in the system. Lowering ep lowers kp which leads to a lower 
density of quasi-particle states at the fermi level which in turn 
reduces the number of mid-gap states in the vortex core. 

Of course, our analysis of the mini-gap in proximity in- 
duced superconductors only accounts for possible quasiparti- 
cle states in the vortex core that were localized at the TI sur- 
face. Therefore the geometry must ensure that there are no 
quasiparticles in the superconductor itself. This can be ac- 
complished by using a discrete vortex as in the tri-junction 
geometry. We believe that the estimates for the continuous 
vortex carry over to the discrete vortex case. 

So far the vortex core size has been taken to be of the order 
of the coherence length in the TI, — — v / A ^ v/X. If 
the vortex core size is taken to be smaller, as is expected for 
a discrete vortex, the numerical calculations lead to an even 
larger AE (Fig. 2). As is clear from Eq. ( |27] l, for [/ < A, the 
excitation gap in a vortex can be of order A, which is A 
for chemical potential on the surface of the TI tuned such that 
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U ^ X (see Eq.fTsTi which can be of the scale of 100 meV (i.e. 
of the scale set by ep). The constraint [/ < A ^ 0.5 meV 
within the Dirac point, that is obtained without the central re- 
sult Eq.([T8]l. is difficult to achieve in experiments because of 
impurity induced density fluctuations as in graphene.^** The 
ability to tune the chemical potential U on the TI surface in- 
dependent of the fermi energy of the superconductor e p leads 
to a significant enhancement over the case of a chiral p-wave 
superconductor where the chemical potential U ^ ep ^ A 
such that the mini-gap from Eq. ( |27l ) scales as AE ^ A'^/ep 
consistent with previous estimates.- 

VII. CONCLUSION 

In conclusion, we have shown that Majorana fermion ex- 
citations in proximity-induced s-wave superconducting sys- 
tems are much more robust to thermal decoherence effects 
than in regular chiral p-wave superconductors. In the latter 



system, the excitation gap protecting the Majorana modes, the 
so-called excitation gap, is given by ~ A'^/ep, which is a 
prohibitively low energy scale ^0.1 mK. On the other hand, 
for proximity-induced s-wave superconducting systems^' 'i, 
which have generated a lot of recent interest^"*"^^, and in the 
case of sufficient transparency of the barriers, the mini-gap 
can be made as high as ~ A ~ 1 K. The possible orders 
of magnitude enhancement of the mini-gap in these systems 
helps bring the observation of non-Abelian statistics to the 
realm of realistic, accessible, temperature regimes in experi- 
ments. Thus, the proposal of Fu-Kane^ and that of Sau et al. — 
with appropriate control of the proximity effect and feature 
sizes appear at this state to provide the most robust platforms 
for the observation of Majorana fermions and the implemen- 
tation of TQC. 
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University start up funds for support. 
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